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1.  INTRODUCTION 


Most  of  the  presently  available  codes  for  computing  the  solution  to  an  electromagnetic  tran¬ 
sient  problem  are  written  to  solve  a  specific  class  of  problems.  If  one  is  interested  in  solutions  to 
several  different  types  of  problems,  one  must  acquire  familiarity  with  a  number  of  different  codes. 
Each  code  must  be  examined  for  its  range  of  validity,  and  it  may  have  to  be  modified  to  fit  specific 
needs.  This  process  may  be  prohibitively  time-consuming,  particularly  if  one  is  interested  in  ap¬ 
proximate  solutions  for  a  wide  variety  of  problems,  rather  than  a  very  accurate  solution  to  a 
specialized  case.  It  is  the  former  type  of  solution,  combined  with  information  about  the  antenna  or 
cable  impedance,  that  often  suffices  for  a  vulnerability  analysis  of  a  given  circuit  attached  to  the 
antenna  or  the  cable. 

The  computer  code  TEMPO  was  developed  with  the  objective  of  providing  approximate  solu¬ 
tions  for  a  wide  class  of  electromagnetic  pulse  (EMP)  interaction  problems  in  a  form  usable  by  a 
circuit  analysis  code.  The  application  of  TEMPO  was  to  be  simple  and  efficient  so  that  the  user 
would  require  a  minimum  amount  of  familiarity  with  the  code,  and  the  code  was  to  be  flexible  so 
that  a  wide  variety  of  solutions  could  be  performed  for  a  number  of  different  conditions.  Finally, 
the  problems  addressed  by  TEMPO  were  to  be  practical  cases  that  could  ultimately  be  related  to 
commonly  occurring  problems  in  an  EMP  vulnerability  assessment. 


2.  CODE  ORGANIZATION 

Since  TEMPO  is  used  in  the  EMP  vulnerability  assessment  of  complex  military  systems,  it  is 
organized  so  mat  its  application  is  an  integral  part  of  such  an  assessment.  To  do  so,  the  EMP 
coupling  source  is  represented  in  terms  of  either  a  Thevenin  or  a  Norton  equivalent  circuit,  and  a 
circuit  analysis  is  performed  on  the  equivalent  circuit  and  its  load.  Since  the  load  circuit  usually 
contains  nonlinear  elements,  the  circuit  analysis  is  performed  in  the  time  domain  by  using 
availab'e  network  analysis  codes  for  circuits  with  nonlinear  elements. 

By  separating  the  EMP  vulnerability  assessment  into  two  parts  consisting  of  the  source 
equivalent  circuit  and  the  load  circuit,  separate  teams  of  investigators  can  be  assigned  to  the  two 
types  of  problems.  Initially,  the  coupling  analysis  can  be  performed  by  using  simple  linear  loads 
for  the  equivalent  circuit,  and  the  load  circuit  can  be  analyzed  independently  from  the  coupling 
source  by  using  simple  driving  sources.  After  the  coupling  source  equivalent  circuit  and  the  load 
circuit  are  adequately  modeled,  they  can  be  combined  for  a  vulnerability  assessment. 

Most  computer  programs  using  a  large  scale  computer  (like  the  Control  Data  Corp.  (CDC) 
6000  series  or  the  IBM  System/360  or  370  family  of  computers)  are  executed  in  a  batch  mode. 
This  mode  requires  that  the  user  prepare  a  number  of  control  cards  for  running  the  program  and 
data  cards  for  supplying  the  input  parameters  to  the  program.  Since  the  process  of  preparing 
such  cards  can  be  time-consuming,  TEMPO  allows  for  an  option  to  generate  the  necessary  infor¬ 
mation  for  running  the  batch  job  using  two  interactive  programs.  These  programs  are  easily  ex¬ 
ecuted  ft  a  remote  teletypewriter  terminal  and  are  designed  to  aid  the  user  to  accurately 
describe  the  problem  that  he  wants  to  solve. 

With  the  interactive  portion  of  initiating  the  batch  job  completed,  the  job  can  be  placed  in  the 
batch  input  queue  at  the  central  processing  site,  and  output  depends  on  the  normal  turnaround 
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time  at  the  computer  center.  By  running  the  programs  in  batch  mode,  there  are  no  limitations  on 
the  size  of  the  programs  that  can  be  executed,  except  for  those  normally  applicable  to  a  batch 
job.  In  principle,  TEMPO  could  initiate  any  pertinent  EMP  interaction  code  previously  run  in  a 
batch  mode  on  a  large  scale  computer. 

The  result  of  the  computation  initiated  by  TEMPO  is  normally  obtainable  in  the  form  of  a  plot 
of  the  transient  response.  The  plot  devices  can  be  a  remote  line  printer,  a  remote  x-y  bed  plotter, 
or  a  remote  teletypewriter.  By  using  the  teletypewriter,  it  is  possible  to  run  TEMPO  and  obtain  an 
output  even  if  the  only  remote  terminal  available  to  the  user  is  a  teletypewriter. 

The  programs  described  in  this  report  were  written  to  run  on  an  IBM  System/370  Model  168 
computer  using  the  time  sharing  option  (TSO)1  for  the  interactive  portions.  They  were  originally 
developed  to  run  on  a  CDC  Model  6700  computer  using  INTERCOM2  for  the  interactive  portions 
Hence,  TEMPO  can  be  run  on  either  system  with  only  a  few  minor  changes  from  one  version  to 
the  other. 


3.  DESCRIPTION  OF  TEMPO  PROGRAMS 

3. 1  Interactive  Programs 

The  present  version  of  TEMPO  is  written  to  provide  approximate  solutions  to  basic  prac¬ 
tical  coupling  problems  for  interactions  with  cables,  antennas,  and  apertures.  There  are  many 
programs  and  subroutines  used  for  performing  the  calculations.  However,  the  user  need  concern 
himself  only  with  the  following  two  interactive  programs  for  initiating  the  computation  process 

3.1.1  Program  INPUT 

Interactive  program  INPUT  allows  the  user  to  select  one  of  five  options.  Three  of  the 
options  (ANTENNA,  CABLE,  and  APERTURE)  lead  to  solutions  of  coupling  problems  by  using 
theoretical  models.  The  CW  (continuous  wave)  option  uses  measured  frequency-domain  data  to 
describe  the  coupling.  The  FIELD  option  does  not  solve  a  coupling  problem,  but  allows  the  user  to 
access  program  DRIVE  without  unwanted  computations.  After  specifying  a  coupling  problem,  the 
user  is  asked  to  supply  information  appropriate  to  that  problem.  For  example,  if  the  choice  is 
ANTENNA,  he  is  asked  to  describe  the  antenna,  specify  the  load  at  the  terminal,  and  add  other 
facts.  The  user  has  an  option  to  obtain  a  brief  synopsis  of  the  model  used  for  the  computation,  in¬ 
cluding  references  for  additional  information  on  the  problem  being  solved.  Once  all  the 
parameters  for  a  given  problem  are  defined,  the  user  is  asked  if  he  wants  output  for  intermediate 
steps  in  the  computation.  These  could  be,  for  example,  the  frequency  spectrum  for  the  transfer 
function,  the  impulse  response,  or  the  driving  field.  Finally,  he  must  specify  the  output  medium  for 
any  intermediate  steps  and  the  transient  response  for  the  problem. 

3.12  Program  DRIVE 

Interactive  program  DRIVE  prompts  the  user  to  enter  the  type  of  incident  held  to  be 
used  for  the  problem.  Several  options  for  the  incident  field  are  available: 


•OS/VS2  TF 0  Command  Language  Reference.  IBM  Technical  News  Letter  GC28  0646  3  ( July  7  9/6) 

1 6000  Series  INTERCOM  3  Reference  Manual.  Control  Data  C orp  Publication  60252800  Rev  B  ( January  797?) 
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a.  A  data  file,  such  as  an  experimentally  determined  incident  field 

b.  One  of  several  types  of  analytic  functions  (a  Gaussian  pulse,  a  sum  of  two  or  four 
exponential  curves,  a  damped  sinusoid,  a  square  wave,  or  a  step  wave) 

c.  A  high  altitude  incident  field  whose  wave  shape  can  have  one  of  the  six  functional 
forms  in  option  b 

d.  An  approximation  to  the  surface  burst  EMP 

Once  the  choice  of  the  incident  field  has  been  made,  there  is  an  option  to  compute  the  driving 
field  either  in  the  time  domain  or  in  the  frequency  domain.  If  the  incident  field  is  determined  from 
experimental  data  in  the  time  domain,  then  there  is  no  choice.  However,  if  the  incident  field  is  an 
analytic  function,  then  computational  error  generally  is  minimized  by  specifying  the  driving  field 
to  be  computed  in  the  frequency  domain. 

3.1 .3  Features  of  Interactive  Programs 

An  attempt  has  been  made  to  anticipate  questions  that  a  user  might  have  during  an 
interactive  session,  and  information  is  provided  to  the  user  in  response  to  a  HELP  command  For 
example,  the  user  may  be  asked  to  specify  the  type  of  coupling  problem.  If  he  issues  a  HElP  com¬ 
mand,  the  program  gives  him  information  on  what  types  of  coupling  problems  can  be  solved  by 
using  TEMPO. 

After  the  interactive  session  is  completed,  a  batch  job  is  generated  This  job  can  be 
submitted  through  TSO  to  the  batch  input  queue  by  issuing  a  SUBMIT' HK3010.  JOB. CNTL'  com¬ 
mand.  Alternatively,  if  a  punched  deck  of  the  job  was  requested  during  the  interactive  session, 
this  punched  deck  can  be  submitted  independently  of  TSO.  Once  a  punched  deck  for  the  |ob  has 
been  obtained,  this  job  can  be  submitted  through  a  card  reader  at  a  batch  terminal. 

3.2  Batch  Programs 

The  remaining  programs  perform  the  numerical  computations  and  generate  an  output 
for  the  problem  specified.  They  consist  of  several  independent  mainline  programs  and 
subroutines  that  pass  pertinent  data  to  one  another  via  temporary  data  sets.  Which  program  is  ex¬ 
ecuted  depends  on  the  type  of  problem  being  solved.  For  example,  the  INVERT  (inverse  Fourier 
transform)  program  may  or  may  not  be  required,  depending  on  whether  or  not  the  computations 
are  performed  in  the  time  domain.  The  following  sections  and  appendix  A  describe  the  mainline 
programs  and  their  subroutines. 

3.2.1  Procedure  10 

Although  not  a  FORTRAN  program,  procedure  10  initiates  all  of  the  files  to  be  used  for 
inputting,  outputting,  and  transferring  data.  This  procedure  creates  a  file  on  unit  2  called  the  10 
file  and  then  stores  the  input-output  control  (10)  data  on  the  file.  The  10  file  is  then  read  by  each 
program.  It  both  passes  pertinent  coupling  parameters  from  one  program  to  the  next  and  controls 
the  task  of  each  program.  The  10  file  also  stores  any  measured  data  submitted  with  the  job 
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3.2.2  Program  FRQ 

Mainline  program  FRQ  reads  .the  10  file,  defines  the  problem  to  be  solved,  and  calls 
the  appropriate  subroutine  to  compute  the  frequency  spectrum  for  the  response.  After  the  fre¬ 
quency  spectrum  is  computed,  the  resulting  effective  height  and  antenna  impedance  are  written 
on  two  separate  files  for  further  analysis.  Subsequent  Fourier  transforms  are  computed  for  2048 
equispaced  frequency  points  by  using  the  fast  Fourier  transform  (FFT)  or  for  256  equispaced  fre¬ 
quency  points  by  using  the  LFILON  transform. 

3.2.3  Program  ZANT 

Program  ZANT  computes  a  lumped  parameter  network  (LPN)  equivalent  to  the  anten¬ 
na  impedance  for  the  problem  being  solved.  The  program  is  used  if  a  time-domain  equivalent  cir¬ 
cuit  is  required  to  represent  the  coupling  problem.  It  applies  when  the  antenna  impedance  shows 
multiple  isolate  resonances,  so  that  the  antenna  impedance  is  easily  represented  as  a  series  of 
parallel  resonant  circuits. 

3.2.4  Program  FIELD 

Program  FIELD  computes  the  incident  electric  field  either  in  the  frequency  domain  or 
in  the  time  domain  as  specified  by  the  user.  The  subsequent  convolution  with  the  impulse 
response  for  the  problem  is  performed  in  the  corresponding  frequency  or  time  domain.  The  op¬ 
tions  available  for  computing  various  types  of  fields  have  been  described  in  the  discussion  of  in¬ 
teractive  program  DRIVE. 

3.2.5  Program  CONVOL 

Program  CONVOL  computes  the  convolution  of  two  functions,  (1)  the  incident  field 
and  (2)  the  impulse  response  for  current  or  voltage  for  a  given  coupling  problem.  For  a  field  scat¬ 
tered  by  an  aperture,  the  second  function  is  the  impulse  response  of  the  aperture.  If  both  func¬ 
tions  are  in  the  frequency  domain,  the  product  remains  in  the  frequency  domain.  If  both  functions 
are  in  the  time  domain,  the  convolution  integral  is  performed  by  using  a  simple  numerical  integra¬ 
tion  in  the  time  domain. 

3.2.6  Program  FOURIER 

Program  FOURIER  transforms  any  one  of  three  types  of  frequency  spectra  from  the 
frequency  domain  to  the  time  domain:  (1 )  the  effective  height  of  the  current  or  voltage  for  a  given 
coupling  problem  or  the  transfer  function  for  the  field  scattered  by  an  aperture,  (2)  the  electric 
field  component  for  the  incident  or  ground-interacted  field,  or  (3)  the  product  of  two  functions 
such  as  the  incident  field  and  the  effective  height  of  an  antenna. 

The  Fourier  transform  is  performed  by  using  either  (1)  theCooley-TukeyJ  algorithm  for 
the  FFT  or  (2)  a  numerical  integration  combining  trapezoidal  and  Filon  techniques  (LFILON).4  The 
FFT  is  used  in  most  cases,  and  the  LFILON  transform  is  available  to  verify  the  results  obtained  by 

*j  W  Cooley  and  J  W  Tukey.  An  Algorithm  tor  the  Machine  Calculation  of  Complex  Fourier  Senes.  Mathematics  of  Computation.  19  ( April 
1965 ) 

*0  G  Dudley.  Numerical  Inversion  of  the  Fourier  Transform  A  Combination  Trapezoidal  and  Filon  JechniQue.  Lawrence  Livermore 
Laboratories.  Livermore.  CA  UCRL  51878  (1975) 
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using  the  F FT.  The  FFT  is  performed  by  using  2048  equispaced  frequency  points,  whereas  the 
LFILON  transform  is  performed  by  using  256  points.  At  present,  the  frequency  points  are 
equispaced  for  the  LFILON  transform.  However,  in  general,  the  LFILON  transform  can  be  used 
with  an  arbitrary  number  of  points  with  arbitrary  spacings  in  the  time  domain,  the  frequency  do¬ 
main,  or  both. 

3.2.7  Program  OUTPUT 

Program  OUTPUT  controls  the  output  processing  for  a  given  run.  The  various 
selected  computed  functions  are  plotted  on  a  bed  plotter,  a  line  printer,  or  a  teletypewriter, 
depending  on  which  device  was  specified  as  an  output  device  The  plot  routines  are  described 
elsewhere  for  the  bed  plotter5  and  for  plots  on  the  printer 6 


4.  METHODS  FOR  EXCHANGING  DATA  BETWEEN  BATCH  PROGRAMS 

Since  TEMPO  is  a  composite  of  various  single  purpose  programs,  the  data  must  be  passed 
from  one  program  to  the  next.  Likewise,  the  instruction  set  and  the  problem  parameters  must  be 
accessed  by  each  program.  The  following  sections  describe  the  input-output  controls,  the 
parameters,  and  the  file  allocations. 


4.1  Definition  of  Temporary  Data  Sets 


Tape  unit 

File  name 

Type  of  data 

1 

— 

Temporary  data  storage  used  by  program  OUTPUT 

2 

10 

Input-output  and  problem  definition  parameters 

3 

FIMP 

Frequency-domain  impulse  response  spectrum 

4 

IMPLS 

Time-domain  impulse  response  spectrum 

5 

FIELD 

Time-domain  electric  field 

6 

PRINT 

Output  to  printer  or  teletypewriter 

7 

SYSOUT=B 

Output  to  bed  plotter 

8 

FFIELD 

Frequency-domain  electric  field 

9 

FRESP 

Frequency-domain  response  to  incident  field 

10 

RESP 

Time-domain  response  to  incident  field 

11 

ZA 

Antenna  and  cable  Impedance 

'T  V  Noon,  Enhanced  Plotting  Software  tor  Use  with  the  Houston  instrument  COMPLOT  Plotter.  Harr y  Diamond  Laboratories  HDL 
TM-75-32  (December  1975) 

• Egon  Marx.  Printer  Version  of  Plots  Made  by  an  Incremental  Plotter.  Harry  Diamond  Laboratories  HDL  TM  75  33  (December  1975) 
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4.2  Description  of  10  Data  Cards 


Card 

Label 

Description 

1 

t 

Time  increment  for  computing  transient  wave  shape 

2 

IOFLAG  (1), 

Flags  to  indicate  which  type  of  output  and  which  type  ot 

1  =  1  to  10 

output  medium  are  desired 

3 

FLAG  (1), 

Flags  to  describe  type  of  problem  to  be  solved  and  type 

1  =  1  to  20 

of  field  to  be  computed,  also  referred  to  as  FLAG  IA. 

IB . IR 

4  to  7 

PARAM  (1), 

Parameters  to  describe  geometry  and  physics  of  scattering 

1  =  1  to  20 

object 

8  to  10 

FPAR  (1), 

1  =  1  to  18 

Parameters  to  describe  incident  electric  field  for  problem 

4  3  Definitions  for  Input-Output  Flags 

Term 

Definition 

MULCURVE 

IOFLG9 

=  1,2,  ,  Depending  on  whether  RESP,  FIELD.  .  .,  contain 

multiple  curves 

NCURVE 

IOFLG10 

=Number  of  curves  -  1  that  are  to  be  plotted  if  MULCURVE  *  0 

FSPEC 

IOFLG4 

=0  plot  magnitude  and  phase 
=1  plot  real  and  imaginary 

IOFLG(l),  1  = 

1.  2.  3,  5,  6,  7,  8.  are  the  outputs  of  files  RESP,  FIELD,  IMPLS. 

FRESP.  FFIELD.  FIMP.  and  ZA.  respectively,  and  are  plotted  in  accor- 

dance  with  I0FLG4  and  the  following: 

Linear  x  and  y  scales: 

IOFlG(l)  =1  plot  output  on  teletypewriter 

=2  plot  output  on  printer 

=3  plot  output  on  plotter 

Allowed  values  for  lOFLG(l)  I  =  1.  2,  3,  5.  6.  7.  8 

Linear  x  scale,  log  y  scale  for  magnitude,  real,  imaginary 
and 

Linear'  x  scale,  linear  y  scale  for  phase: 

IOFLG(l)  =5  plot  output  on  printer 

=6  plot  output  on  plotter 

Allowed  values  for  IOFLG(l)  I  =  5,  6.  7,  8 

Log  x  scale,  log  y  scale  for  magnitude,  real,  imaginary 
and 

Log  x  scaie,  linear  y  scale  for  phase 
IOFLG(l)  =8  plot  on  printer 
=9  plot  on  plotter 

Allowed  values  for  IOFLG(l)  I  =  5,  6.  7.  8 
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4  4  Flag  Card 


The  Uag  card  has  20  integers  in  columns  1  to  20:  flags  IA  to  'T.  'rhe  flags  control  ihe 
routing  within  the  mainline  programs  and  the  subroutines  so  that  the  desired  problem  will  be  solv¬ 
ed.  Figures  1  to  8  indicate  the  uses  of  the  flags. 


Figure  1.  Definition  of  flags  1C,  ID  (mainline — INPUT). 


Figuie  2. 


Definition  of  flags  IA,  IB,  1C,  IE,  IF  (CABLE) 


2 


Figure  3.  Definition  of  flags  IA,  IE,  IF,  IN  (ANTENNA) 
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Figure  6.  Definition  of  flags  IK,  IM.  IN  (mainline— DRIVE) 
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Figure  8  Definition  of  flag  IL  (SIMULATOR) 


The  following  flags  are  used  in  addifon  to  those  defined  in  figures  1  to  8 
Flag  Description 

IH  Computation  of  equivalent  network  for  antenna  impedance 
IH  *  0  computation  not  required 

IH  *  1  computation  required 

IS  Method  used  for  computing  Fourier  transform 

IS  =  0  use  FFT 

IS  =  1  use  LFILON  transform 

IT  Debug  option 

IT  *  0  no  debug  option 
IT  *  1  debug  option  desired 

4  5  Definition  of  ARRAY  Parameters  (PARAM)  1  to  20 
Parameter  Definition 


PARAM(  1 1  «  RHO  (resistivity)  definition  when  solving  cable  problem 

RHO  =  -  5E-5  if  analytical  model  Cl  is  used 
RHO  *  1  /SIGMA  (SIGMA  =  ground  conductivity)  if  analytical 
model  C2  is  used 


PARAM(2) 


PARAMO) 


PARAM(4) 


PARAM(5> 


For  cable  problem.  =  R1  of  senes  resistance-inductance  capacitance  (RLC) 
at  cable  terminals  (ohms) 

For  antenna  problem,  *  R  of  RLC  load  at  antenna  feed  point  (ohms) 

For  continuous  wave  (cw)  problem.  =  R  *  load  resistance  (ohms) 

For  aperture  problem.  =  H  =  length  of  aperture  (meters) 

For  cable  problem.  =  LI  of  series  RLC  at  cable  terminals  (henries) 

For  antenna  problem  =  L  of  RLC  load  at  antenna  feed  point  (henries) 

For  cw  problem.  =  IND=  load  inductance  (henries) 

For  aperture  problem.  =  W  =  width  of  aperture  (meters) 

For  cable  problem,  =  Cl  of  senes  RLC  at  cable  terminals  (farads) 

For  antenna  problem,  =  C  of  RLC  load  at  antenna  feed  point  (farads) 

For  cw  problem,  »  C  =  load  capacitance  (farads) 

For  aperture  problem,  =  R  =  distance  from  aperture  for  computing  electric 
field  (meters) 

For  cable  problem  =  R2  of  series  RLC  at  cable  terminals  (ohms) 

For  antenna  problem, 

*  H  -  length  of  antenna  for  monopole  and  dipole  antennas  (meters) 

*  B  =  loop  radius  fci  loop  antennas  (meters) 

»  H  =*  missile  length  for  missile  problems  (meters) 

=  Hi  *  0  3  for  AS1852  antenna.  *  0  35  for  AS1852  antenna  (meiers) 
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Parametei 


Definition 


PARAM(6) 


PARAM(7) 


PARAM(8) 


PARAM(9) 


PARAM(IO) 


PARAM(II) 


PARAM(12) 


For  cable  problem,  =  L2  of  series  RLC  at  cable  terminals  (henries) 
For  antenna  problem, 

=  D  =  antenna  diameter  for  dipoles  or  monopoles  (millimeters). 

=  R  =  wire  radius  for  loop  antennas  (meters). 

=  A  =  missile  radius  for  missile  problem  (meters). 

For  cw  problem,  =  RES  =  antenna  resistance  (ohms). 


For  cable  problem,  =  C2  of  series  RLC  at  cable  terminals  (farads). 

For  antenna  problem, 

=  X  =  H*XIN  (H  =  missile  length)  representing  relative  position  between 
missile  tail  and  missile  nose  for  current  calculations  (XIN  =  decimal 
number  between  0.  and  1 .). 

=  X  =  HP»XIN  (HP  =  missile  plume  length)  representing  relative  position 
between  plume  tip  and  missile  tail  for  current  calculations  (XIN  =  decimal 
number  between  -1.  and  0.). 

=  LL  =  2000  for  system  antennas  AS1852  and  AS1853 

For  cable  problem,  =  L  =  cable  length  (meters). 

For  antenna  problem, 

=  HP  =  plume  length  for  missile  problem  (meters). 

=  ZOI  =  50  for  system  antennas. 

For  cable  problem,  =  X  =  value  between  0.  and  1.  representing  fractional 
distance  along  cable  for  computing  response. 

For  antenna  problem, 

=  R1  of  series  RLC  termination  on  cable  connected  to  antenna  for  system 
antennas  (ohms). 

-  SIGMA  =  plume  conductivity  for  missile  problems. 

=  1.E7  (automatically)  if  plume  length  =  0. 

For  cable  problem,  =  H  =  height  of  cable  aboveground  (meters). 

For  antenna  problem, 

=  HG  =  height  of  antenna  terminal  aboveground  for  system,  loop,  monopole, 
and  dipole  antennas  (meters) 

=  HG  =  height  of  missile  aboveground  for  missile  problem  (meters). 

For  cw  problem,  =  HG  =  height  of  test  point  aboveground  (if  antenna 
resistance  is  greater  than  zero)  (meters) 

For  aperture  problem,  =  HG  =  height  of  center  of  aperture  aboveground 
(meters). 

For  cable  problem,  =  D  =  conductor  diameter  (meters) 

For  antenna  problem, 

=  LI  of  series  RLC  termination  on  cable  connected  to  antenna  for  system 
■  antennas  (henries). 

=  ALPHA  =  damping  constant  for  missile  problems. 


For  cable  problem,  =  S  =  separation  distance  between  cables  for  two- 
conductor  cable  (meters). 

For  antenna  problem, 

=  RP1  -  inner  plume  radius  for  missile  problem  (meters). 

=  Cl  of  series  RLC  termination  on  cable  connected  to  antenna  for  system 
antennas  (farads). 
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Parameter 
PARAM(1 3) 

PARAMO  4) 

PARAMO  5) 

PARAMO  6) 

PARAMO  7) 

PARAMO  8) 

PARAMO  9) 
PARAM(20) 

4.6  Definitions 


Definition 


For  antenna  problem,  =  outer  plume  radius  for  missile  problem  (meters). 
For  frequency  plot,  =  multiple  of  frequency  increment  (DF). 

For  frequency  plot,  *  maximum  frequency  for  response  plot  (hertz). 

For  frequency  plot,  =  maximum  frequency  for  field  plot  (hertz). 

For  frequency  plot,  =  maximum  frequency  lor  impulse  plot  (hertz). 

For  frequency  plot,  =  maximum  frequency  for  impedance  plot  (hertz). 

For  time-domain  output,  =  maximum  time  for  transient  output  (hertz). 

For  time-domain  output,  =  maximum  time  for  field  plot  (hertz). 

For  time-domain  output,  =  maximum  time  for  impulse  plot  (hertz). 

of  Field  Parameters  (FPAR)  1  to  18 


Parameter  Definition 


FPAR  (1)  For  function-described  field,  =  peak  amplitude  (volts/meter)  =  VMAX. 

For  simulator  field,  =  YMAX  =  peak  free  field  at  50  m. 

For  surface  field,  =  YMAX. 

FPAR  (2)  For  Gaussian  field,  *  WIDTH  *  pulse  width  (nanoseconds). 

For  two-exponential  field,  =  ALPHA  =  one  of  two  decay  constants  (per 
second). 

For  sine-function  field,  =  OMEGA  =  (2*PI)/PERIOD  *1.E9  (PERIOD  in 
nanoseconds). 

For  four-exponential  field,  =  ALPHA  =  one  of  four  decay  constants. 

For  square-wave  field,  =  WIDTH  =  pulse  width  (nanoseconds). 

For  surface  field,  =  ALPHA. 

FPAR  (3)  For  two-exponential  field,  =  BETA  =  one  of  two  decay  constants  (per  second). 
For  sine-function  field,  *  ALPHA  =  1  ,E9/e-folding  time  (e-folding  time  in 
nanoseconds). 

For  four-exponential  field,  =  BETA  =  one  of  four  decay  constants. 

For  surface  field,  =  BETA. 

FPAR  (4)  For  four-exponential  and  surface  field,  =  GAMMA  =  one  of  four  decay 
constants. 

FPAR  (5)  For  four-exponential  and  surface  field,  =  ETA  =  one  of  four  decay  constants 

FPAR  (6)  For  fou'  exponential  and  surface  field,  =  Al  =  one  of  the  three  "weighting" 

constants. 

FPAR  (7)  For  four-exponential  and  surface  field,  =  A2  =  one  of  three  "weighting" 
constants. 
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Parameter  Definition 

FPAR  (8)  For  four-exponential  and  surface  field,  =  A3  =  one  of  three  "weighting'' 
constants 

FPAR  (9)  For  simulator  field,  =  RO  =  radial  distance  from  simulator  (meters). 

For  high  altitude  field,  =  THETA  -  angle  of  incident  field  with  respect  to 
ground  (degrees). 

FPAR  (10)  For  simulator  field,  =  R2  =  height  of  simulator  aboveground. 

For  all  other  fields,  =  PHI  =  polarization  angle  of  total  E-field  vector  measured 
up  from  horizontal  (degrees). 

FPAR  (11)  For  computing  reflected  component  of  high  altitude  incident  field,  -  SIGMA 
=  ground  conductivity  (default  value  =  0.01  mho). 

FPAR  (12)  For  computing  reflected  component  of  high  altitude  incident  field, 

=  EPSILON  =■=  ground  dielectric  constant  (default  value  =  15). 

FPAR  (13)  For  following  three  cases  (loop  antennas  only),  FPAR  (13  to  18)  have  following 
values. 

For  loop  horizontally  oriented  with  THETA  (angle  of  field  with  respect  to 

ground)  less  than  or  equal  to  45  deg  (orientation  of  E-field's  polarization  is 

unimportant  here). 

For  loop  vertically  oriented  with  E-field  vertically  polarized  parallel  to  plane  of 

loop  (THETA  must  be  less  than  or  equal  to  45  deg). 

For  loop  vertically  oriented  with  E-field  horizontally  polarized  (no  restrictions 
on  field's  orientation  with  respect  to  loop  (THETA  must  be  greater  than  45 
deg). 

FPAR  (13)  =  EX  -  0.  FPAR  (16)  —  L  =  — 1. 

FPAR  (14)  =  EY  =  1.  FPAR  (17)  =  M  =  0. 

FPAR  (15)  =  EZ  =  0  FPAR  (18)  =  N  =  0. 

For  following  three  cases  (for  loop  antennas  only),  FPAR  (13  to  18)  have 
following  values: 

For  loop  horizontally  oriented  with  THETA  greater  than  45  deg  (orientation  of 
E-field's  polarization  is  unimportant). 

For  loop  vertically  oriented  with  E-field  vertically  polarized  perpendicular  to 
plane  of  loop  (THETA  must  be  greater  than  45  deg) 

For  loop  vertically  oriented  with  E-field  horizontally  polarized  (no  restrictions 
on  field's  orientation  with  respect  to  loop)  (THETA  must  be  less  than  or  equal 
to  45  deg) 


FPAR  (13)  =  EX  =  0 
FPAR  (14)  «  EY  =  1 
FPAR  (15)  =  EZ  =  0 


FPAR  (16)  =  L  «  0. 
FPAR  (17)  =  M  =  0. 
FPAR  (18)  =  N  =  -1 
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APPENDIX  A.— PROCEDURES  FOR  COMPUTER  CODE  TEMPO 
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APPENDIX  A 


A-1.  PROCEDURE  TEMPO 

Procedure  TEMPO  (listing  A-1)  can  be  used  through  the  time  sharing  option  (TSO)  to  execute  in 
teractive  programs  INPUT  and  DRIVE.  After  procedure  TEMPO  is  completed,  the  new  batch  job  can  be 
submitted  to  the  job  stream  by  specifying  a  SUBMIT  JOB  command.  The  various  data  sets  such  as  10. 
CW,  and  JOB  can  be  examined  and  altered  by  using  the  TSO  EDIT  command 

Listing  A-1.  Procedure  TEMPO 


00010  P ROC  0 
00020  CONTROL  NOMSQ 

0  0030  FRE  E  FILE(FT01FOQ1,FTQ2FOQ1  , FT03F00  1  ,  FTQ4F00 ] ,FT05F001>FT07FQQ1) 
00040  CONTROL  MSG 

00050  ALLOCATE  OS  NAME (  *)  FILEC  FTO 1 F0  0  1  ) 

00060  ALLOCATE  DSNAMEC*)  F I LE ( FTO 2F0 0 1 ) 

00065  ALLOCATE  DSNAMECIO)  F I LE< FTO 3F0 0 1 ) 

00070  ALLOCATE  DS  NAMEC  JOB  .CNTL)  F ILEC FT04F0 0 1 ) 

00080  ALLOCATE  DSNAMEC  CWDAT)  F I LE  C  FTO  5F0  0 1 ) 

00090  ALLOCATE  SYSOUTCB)  FILEC FTO 7F0 0 t ) 

00100  FORTGO  ( I NPT ) 

00110  FORTGO  (DRIVE) 

00120  CONTROL  NOMSG 

0  0  130  FREE  FILEC  FTO 1 F00 1 , FTQ2F0  0 1 > FTO 3F0 0 1 , FTO 4F0  0 1 , FTO  5F0 0 1 , FT07F0  0 1 ) 
00140  CONTROL  MSG 
END  OF  DATA 
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A-2.  PROCEDURE  IO 

Procedure  10  (listing  A-2)  is  used  the  firs!  time  that  the  10  parameters  are  to  be  read  from  cards  A 
the  same  time  aa  temporary  data  sets  to  be  used  in  a  TFMPO  batch  |ob  are  defined 

Listing  A-2  Procedure  ... 


//13H  PRLC 

//*  T3R  HARRIS,  303  — b  99  —  1  7  4  N  ,  1  JU  l  Y  1  ?  8  d  ,  INDEFINITE 
//*  IV1T1ATE5  INPJT  AND  DATA  FIEFS 
//» 

/  /  0 1 H I  S  PRE'CEJJRE  READS  IN  10  DATA  CAROS  USED  TU  Of  SCK  1  BE  THE  PRUblFM 
//#  TU  BE  EXEtJTED  eY  TEMPO.  THIS  INCLUDES  DT.IUFLAGU),  FLAGll), 
//»  PARAil  1  )  .  AND  FPAR  ID 
/!  o 

II C3PY  EXtC  P  C  M  = IEHGENEK 
II SY51N  DD  DUMMY 
//SYSPR1NT  Dt  DU  MR  Y 

//SYSJT2  DD  D SN^EClotDISP-I NEW, PASS), JNIT=SYSDA, 

//  SFA(F-(TRK,I1,11), 

//  DCBMRECFM  =  FH,lKKl--feO,EttKSin--8  0C) 

//FT(j3fDCl  DO  DSN-U.F  IMP, DISF  MNEw, PASS). 

II  UNIT =V  13, SPACE  =  (TRK,  (1  ,  1  )  I, 

//  CCbMRECFM  =  VS,BEKSTZf*4(,9fc) 

//FTUAFJOl  DD  DSN-LCIMPL  S,D1 SP  -  <NEW,PASS)  , 

//  li  N I  I  -  V  13, SPACE  -  (TRK,(1  ,1  )  ), 

II  DCb  =  I  RECF  H=VS,BLKSIZE  =  4  09  6) 

//FI0SFD01  DD  DSN  =  GLF1F  LD.DI SP  -  (NEW, PASS)  , 

//  UNIT-V1U.SFACE  =  I  TRK  ,  (  1  ,  1  )  ), 

//  DCb  =  <  RECFM  =  vS.BLKSUE=<.r961 

//FT06F301  CD  OSN-ttPRINT,DlSP-(NEW,PASS)  , 

//  UN1T  =  Y13, SPACE  =  TTRK,  19, c)  ), 

//  DCb-IRECFMrFA,LRECE-]33,PLKSWt  -133) 

//FT07F001  DD  DSN=LCPUNDATA»DISP=  IN: W,P  ASS), 

II  JMT=V1L, SPACE  =  HRK  ,11  ,1  )), 

//  OCB-(RECFM  =  FB,LR  tCL  -8U.BLKS1  Zl  =  b  30  I 

//FT03F301  DO  DSN-ttFF  1ELD,D1SP  =  t  NE  *  «  P  A  S  S  1, 

//  UNI I=VI3. SPACE  =  (TKK,C1,I  I  ), 

//  DCb-IRECFM=VS,6lKSI2E  =4090 

//Fr09F301  DC  DSN  =  ttFPK[)D,DISP  -  (NEW, PASS)  , 

//  UNII  V  13, SPACE  =  t  TKK  ,  (1,1  )  ), 

//  DC B  =  ( R E C E M - V 5  ,  B L KS  1  Z E = 4  9b) 

//FT13FD01  DJ  DSN  =  ILRE SP  ,D ISP  =  INfcW, ■> ASS >  , 

//  U NIT  =  V1D, SPACE  =  ITRK, (1,1)), 

//  0Cl-IRECFM-VS,BLKSIZF-4(9fr) 

//FT  1 1 F  301  DD  OSN=CCZA , D 1 S  P  = <NE W . PASS ) , 

//  UNIT=V13,SPACf  MTRK.U,!)  ), 

//  CCE=(KECFM-VS,BIKS1ZE=4  96) 

// FT20F001  DD  DSN-tCNL JN,D)SP=INEW,PASS  ) , 

//  UN  1 1 = V ID,SPAC£-(TRK, 19,9)  ), 

//  DCP-«AECFM=WS,BLKS1ZF=4( 9t) 
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A-3.  PROCEDURE  101 

Procedure  10 1  (listing  A-3)  reads  cards  into  a  specified  temporary  data  set  The  default  value  for 
the  temporary  data  set  is  10,  since  this  procedure  is  normally  used  to  read  a  new  set  of  parameters  into 
the  temporary  data  set  10 

Listing  A-3  Procedure  101 


//iai  PR  3  C  EILE-'ltlO'.CUT-X 

//»  M  .  J.  STARK,  664-6234  ,  17  AJG  77,  INDEFINITE 

//*  READS  13  DATA 

/  /  1  0  I  tXEC  PGM=  IE  f-GE  N£  R 

//SYSIN  DO  DUMMY 

//SYSPTINT  DU  DUMMY 

//SYSJT2  DC  OSN-tE ILE , D I SP  =  (GLD,PASS) 

/ / C  3PY  EXEC  PGM  = IE  fcGENEft  .COND-EVE  A 
//SYSIN  UD  DUMMY 
//SYSPRINT  DL  DUMMY 

//SYSJT1  DC  L>SN  =  £.F  1LE  ,D1SP  =  I0LC,PASS) 

//SYSJT2  DO  SYSllUT  ^EDUT 

A-4.  PROCEDURE  102 

Procedure  102  (listing  A-4)  copies  a  short  term  or  permanent  data  set  into  a  specified  temporary 
data  set  The  name  of  the  input  data  sef  must  be  specified  on  the  EXEC  statement.  The  default  value  for 
the  temporary  data  set  is  FIMP.  since  this  procedure  is  normally  used  to  copy  continuous  wave  (cw)  data 
into  the  temporary  data  set  FIMP 

Listing  A-4  Procedure  102 

/  /  1 3  2  PR  DC  1N= ,UJT - *IEF IMP* 

//•  *.  J.  STARK,  664-6234,  2  APR  77,  INDEFINITE 

//*  COPIES  FILES 

/ / C  3  PY  EXEC  PGM-IEfcGENFR 

//SYSIN  DD  DUMMY 

//SY5PRINT  DO  DUMMY 

//SYSJU  DD  DSN=tl  N.DISPMOID.KEEF  I 

//SYSJT2  DO  USN-tauT.OISPMCiLD.PASS) 

//  EXEC  PGM= IE  bGENER  ,COND=EVEN 

//SYSIN  DD  DUMMY 
//SYSPRINT  DD  DUMMY 

//SYSUT1  DD  DSN=COUT ,DISP - (OLD, PASS  I 
//SYS- T  2  DD  SYSOUT=A 
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A  S.  PROCEDURE  FRQ 

Procedure  FRQ  (listing  A-5)  loads  and  executes  program  FRQ  witn  the  appropriate  subroutine  The 
appropriate  FAMP  subroutine  must  be  specified  on  the  EXEC  statement  Otherwise,  the  default  value  of 
FAMP  1 234  is  used  It  more  than  one  frequency  spectrum  is  to  be  computed,  then  FDISP  =  MOD  must  be 
specified  on  the  EXEC  statement 

Listing  A-5  Procedure  FRQ 

//FREQH  PRLC  FAMP=FAMP123A,F01SPrLLb,EDlSP=  ULD 

//*  TDM  HARR  1 S, 303-599-1  745,  1  JULY  1x80  ,  INDEFINITE 

//»  CDMPt'T  ES  FAMP  FILES 

//IKED  EXEC  PCM  =  IE  ML  ,RE  0  IIJN  =  2  50  K  .  PAR  M  =  'NCMAP  ,U  ST  .NLXftEF  .LET  * 

//« 

//*  THIS  PROCEDURE  EXECUTES  THE  FRrOJENCY  IMPULSE  RESPUNSE  MAINLINE 
//*  PRJCRAM.  THE  PROCEDURE  CARD  MJST  HAVE  FAMP^F AMPSUBROUT INE  TU  BE 
//»  EXECJTED. 

//* 

//Sysprint  dl  Dummy 

//STSLI3  DO  0  SN= 'HK301A.TEMPDA  * .DIS^-SHR 

//  DO  DSN=SYS1.F0RTL IB,D1SP=S7R 

//  DO  D  S  N  =  SYS1 .PACLOAO  «  0 I SP  =  S  H  R 

//SYSJT1  DO  UNIT =SYSDA, SPACE M IC24.I  200  ,201) 

//SYSLMOO  DO  DSN=CU.0SETlMAINl,UMT=5rSDA,DISP=«, PASSE, 

//  SPACE:(TRKt(10,10,2ll 

//SYS11N  DO  DSN-HK301A .LKED1N4I CF  AM  J I , D ISP'S HR 
//FRED  EXEC  PGM-*.LKED.SYSLMtlD,CONDz  E8,LT  ,LKEDI 
//FT02FD01  DD  0SN  =  1 1 1  (1 ,  D  1  SP  -  « Ul  D  .  PA  S  S  I 
//FT03FD01  DD  DS  N  =  l CF IMP  ,D  I  S  P  - ( CFD1SP,PASS1 
//FT  ObF  001  CD  DSN  =  C CPU  I  NT , D I SP - I  MOO , P A S S 1 
//  F  I  08  F  00  1  DO  DSNsCCF  F  1 1  LD  ,D  1  SP  MIX  J  ,PASS  I 
//FT09P001  DD  DSN=CCFPkODtDISP=»OLD,PASS> 

//F  T  1 1  F  001  DD  USN  =  EUA,DISP=  » t  Z  D  I  SP  ,  P  AS  S  » 

//FI2DF001  DD  DSN = t CNL 1 N  ,D 1 SP  =  t LL D, P ASS > 

//DD 1  DD  DSN-llOOSET ,DISP=lPlb, DELETE) 
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A-6.  PROCEDURE  EFIELD 

Procedure  EFIELD  ('isting  A-6)  loads  and  executes  programs  EFIELD  and  EFIEL1  If  more  than  one 
field  is  to  be  computed.  FRQ  =  MOD  or  TIM  =  MOD  must  be  specified  on  the  EXEC  statement,  depending 
on  whether  the  field  computation  is  in  the  frequency  domain  or  the  time  domain 

Listing  A-6.  Procedure  EFIELD 

//EFIELOH  Fftt'C  F  R3  =OLD , T  1  M  =L1  L  D 

II*  TQM  HARRIS, 303-599-1 74S,  1  JULY  1J8D  ,  INDEFINITE 

//*  COMPUTES  EFIELD 

II* 

//*  THIS  PAQC  EDURE  EXECUTES  THE  EFIELD  MAINLINE  PROGRAM. 

II* 

/ / L  K  ED  EXEC  PGM  = I E ML ,R E G I  ON  =  3S0K . PA  * M = • NDMAP .1  I  ST ,NDXRE F ,L E T • 

//SY SPRINT  00  DUMMY 

//SYSLI3  DO  0SN=  ‘HK3014  .TEMPCA »,D1S?=SHR 
//  DO  0SN=SYS1 .FLRTL lb ,  D  I SP-SYR 

//  DO  DSN=SYS1.PACLUA0,0!SP=SyR 

//SYSJT1  00  JNIT=SYSOA, SPACE-11'  2A,I3.0  ,?n)I 
//SYSLMQD  DO  DSN=CIGUSE  T  (MAIN  I  ,UM  US  YSD  A ,  D1  SP=  I  .PASS  >  , 

//  SPACE-(TRK,(10,10,2)I 

//SYSLIN  DC  DSN=HK301'..LKE0IN4(EF  IE  ,  D  I  ,  DI  SP  =  SHR 
// E  F  LD  EXEC  PGM=*.LKEO.SYSLMCiD,COND  -  (B  ,LT  ,LKFD  1 
/  /  F  T  03  F  DO  1  00  DSN-tClU,DlSP-(ClLO,PASSI 

//FT03FD01  DO  DS  N  =  CIF IMP, DISPr(OLC, PASS  t 
//FTOSFDOl  CO  DSN  =  l IF IE  LD,Dl SP=(ITIM»PASS  1 
// F  T  Ob  F  00 1  DO  OSN^ELPR 1NT.D1 SP * (MU)  , PAS S ) 

//FTOSFDOl  DO  DSN-tlFF ltLL,DISP=(EFTC,PASSI 
//0D1  DD  DSN  =  CCS0SET,01SP=(LltD,CELEIEl 
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A-7.  PROCEDURE  INVERT 

Procedure  INVERT  (listing  A-7)  loads  and  executes  program  INVERT  and  the  appropriate 
subroutine  The  default  value  for  this  subroutine  is  FFT1.  Otherwise,  FOURIER  =  LFILON  must  be 
specified  on  the  EXEC  statement  If  the  Impulse  response,  the  field,  or  the  transient  response  is  to  be  com¬ 
puted  more  than  once,  then  IDISP.  FDISP,  or  RDISP,  respectively,  must  be  set  equal  to  MOD  on  the  EXEC 
statement. 

Listing  A-7.  Procedure  INVERT 

//INVERTH  PRL't  lDlSP*OLD,FDISP=DLD,RDlSP*OLO,FaRlER=FFTl 
//»  TOM  HARRIS, 303-S99-1745,  1  JUIY  1983  ,  INDEFINITE 

//*  COMPUTES  EFT 
//* 

//«  THIS  PROCEDURE  EXtCUTES  THE  F0JR1ER  TRANSFORM  MAINLINE  PROGRAM. 

//*  AS  OF  1  JULY  ,198  0,  THIS  COMPUTES  AN  INVERSE  FOURIER  TRANSFORM, 

//*  FREQUENCY  TO  TIME  DOMAIN,  ANt  IT  IS  EXPECTED  TO  COMPUTE  A  FORWARD 

//»  F0JR1ER  TRANSFORM  TIME  TO  FREQUENCY  DUMAIN,  IN  THE  NEAR  FUTURE. 

//* 

//IKED  EXEC  PGM= IEWL,REG1lN=?S0K,PARM- 'NO MAP, LIST ,NUXREF,LET* 

//SVSPRINT  DO  DUMMY 

//SYSLIB  DO  DSN-  *HK3014.TEMPtl4  *  ,  D  IS  5  =  SH  R 
//  Dl  DSN=SYS1.FURTL IB.DISP-SHR 

//SYSUTl  DO  UNIT  =SYSDA, SPACE -U  24,(20, 2011 

//SYSLM3D  DO  DSN=£.LGOSET  IMA  IN  )  ,UM  TrSYSD  A,D1SP=I  .PASS!  , 

//  SPACEMTRK, <10, 10,2)1 

//SYSLIV  DO  DSN  =  HK3014.LKFD1N4|£.FCAIE»)  ,D1SP=SHR 

//FOUR  EXEC  PGM-«.LKED.SYSLM0D»10N0= 1B,LT  ,  LK  ED  1 
//FT0/F001  DO  0SN=CLI0,0ISP=<CILP«FAS5) 

//FT03FD01  CD  DS N -E CF I MP ,D I S F  =  I  OLD, » A SS ) 

//FTORFDOl  DD  DSN-LtIMPLS,DISP-(C  IDISP, PASS) 

//FTuSFDOl  DD  DSN  =  UF  I  E  LD ,  D I  SP  -  UFD  I  SP  ,  P  A  SS  I 

/  / F I  Ob  F  001  DD  DSN-LCPRINT,0ISP  =  1M0D,PASSI 

//FT08FD01  DD  DSN-LLFE IELO,DISP  =  lULD,PASS  ) 

//FT09F001  DD  DS N EPR OD ,D 1 SP  =  I  OLD , P A S SI 

//Ft  10F0OI  DD  DSN  =  tLRESP,DISP-l ERDISP.PASS) 

//DD 1  OD  DSN=tCGOSETfDISP-IOLD,OFLEfEI 
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A-8.  PROCEDURE  CONVOL 

Procedure  CONVOL  (listing  A-8)  loads  and  executes  program  CONVOL  It  more  than  one  convolu¬ 
tion  is  to  be  performed,  then  FRQ  =  MOD  or  TIM  =  MOD  must  be  specified  on  the  EXtC  statement,  depen 
ding  on  whether  the  convolution  is  in  the  frequency  domain  or  the  time  domain 

Listing  A-8  Procedure  CONVOL 

//C3NV3LH  PRJC  T1M-LLD»FRC-ULD 

//*  TDM  HARR  I  St  3  03—599  —  1  7  A  *> ,  1  JUIY  1  J  8  3  ,  INDEFINITE 

//*  C3MPUT  f  S  CONVOLUTION 
//* 

//*  THIS  PROCEDURE  EXECUTES  THE  CQN V3C  UT  ION  MAINLINE  PROGRAM. 

//®  IT  C3MPJTES  EITHER  A  TINE  DOMAIN  C  ON  VOL  UT 1  UN  UF  THE  IMPULSE  RESPONSE 

//*  AND  EFIELD  DRIVE,  LiR  A  FREQUENCY  DOMAIN  PRODUCT  OF  THE  FREQUENCY 

//«  OUNAIN  RESPONSE  AND  FREQUENCY  DOMAIN  EFIELD  DRIVE. 

//* 

/  / 1  R  E  D  EXEC  P  G  M  =  I  ERE  ,H  E  u  ItlN-DSOF  .  PA  R  M  =  •  N  0  M  A  P  ,  L  I  ST  ,NUXREF  ,LE  T  * 

//SYSPRINT  OP  DUMMY 

//SYSL13  DO  0SN=  *HR301A  .Tl  MP(:A  •  ,D  ISP  =SHR 
✓  /  DO  9SN=SYS1  .F  PR  TL  lb  ,  D  I  SP-SiR 

//SYSJT1  DD  JN IT =SYSDA  .SPACE =1  li  2A,I 2  DO  ,20) I 
//SYSLM3D  DD  DSN=  E.E.GC1SE  T  (MAIN  I  ,  UN  IT  =  S  YSD  A,  D  I  SPM  ,P  ASS  I  , 

//  SPACE = ITRK  ,  (1 u,  IP  ,i  I  I 

//SYSL1N  DD  DSN=HK301A.lKf  Dl N4IC0NV31 I , D I SP  =  SHR 

/ / C  3NV  EXEC  PGM=«.LKED.SYSLMUD,COND= EB.LT  ,LKFD) 

//FT02F001  DD  DSN^LCIO , D 1 SP  = (UL  C  .PASS) 

//rr03F0Ul  DD  DSN  =  ttF I KP  ,D I S P = ILL D, P ASS  I 
/ / f T0AF3D1  CD  DSN-tLIMPlS,Dl SP  s  C  OLD  ,  P A  S  S  I 
//FIOSFOOl  DD  0SN-E.EPIELD,DISP-(01D,PASSI 
//FI  Ob  F  301  DD  DSN-CCPR INT ,DI SP  * I  M  GO  ,  P  A  S  SI 
//FT06F001  DD  OSN-ttFF IE ID.OISP - EU.D,PASS  I 
//FT09F301  DD  DSN-ICFPROD,OI SP  -  ILFRj,PASS  I 
//FT  1 0  F  001  DD  DSN  =  ttRESP,DISP  =  UT  IM.PASSI 
/  /  DD  1  DD  DSN^LLGOSE T ,D 1 SP= IDLD .DELE T E I 
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A-9.  PROCEDURE  OUT 

Procedure  OUT  (listing  A-9)  loads  and  executes  program  OUT  The  default  value  ot  O-FPLOT  = 
HKZ001  OFFPLOT'  allows  the  user  to  write  output  to  the  bed  plotter  on  magnetic  tape  for  subsequent  off 
line  plotting  The  substitutable  parameter  OFFPLOT  must  be  set  equal  to  SYS1  DUMMYLIB  for  on-line 
plotting  The  default  value  of  HLDPLT  =  YES  allows  the  ou>  '  ,o  the  plotter  to  be  placed  in  hold  so  that  the 
operator  can  mount  a  magnetic  tape  for  the  off-line  plot  option  Otherwise,  the  user  must  specify  HLDPLT 
=  NO  on  the  EXEC  statement  Substitutable  parameters  SYST  and  HLDPFtT  are  used  when  a  job  is  submit¬ 
ted  through  TSO  and  it  is  desired  to  recover  the  output  at  the  TSO  terminal  such  as  a  plot  on  the 
teletypewriter  In  that  case,  the  user  must  set  SYST  =  2  and  HLDPFtT  =  YES  The  output  can  then  be 
recovered  at  the  TSO  terminal  by  issuing  the  OUTPUT  (JOBID)  command. 

Listing  A-9  Procedure  OUT 

//DOTH  PRUC  OFFPLOT- •  H  K  30  1 4  .OFFPLOT  *,SYST-A,HLDPLT-YES,HLDPRT=NO, 

//  01 S-PASS.PUN-B 

//«  TOM  HAPR  1  S,  3U3-S99-1  74S,  1  JULY  |98  j  ,  INDEFINITE 

//•  OJTPlfT  RL  Jf  INE 

//• 

//•  THIS  PROCEDURE  PLOTS  THE  UUTPJI  RESULTS  FROM  TEMPO  AS  CONTROLLED 

//•  B Y  THE  I0FLAGS  IN  THE  10  DATA  DECK 

//• 

//LKEO  EXEC  PGM- IE WL,KEGK,N  =  2S0K .PARK- 'NCMAP ,L 1ST .NOXREF .LET • 
//SYSPRINT  01‘  DUMMY 

//SYSL18  DC  DSN-  CUFFPLO  T .DISP-SHR 
//  DO  DSN  *0X3010. ANAPAC * .OISP-SHR 

//  00  DSN  -  SYSl  .F  URTL  IE>  ,0  I  SP  -SHR 

//  OD  DSN-SYS1  .PAGLL1AD  , D  I  SP  =S  -tR 

//  UD  DSN^'HXJUIA.TEMPGA'.DISP^SHA 
//SYSJT1  01  JN  IT  -SYSDA, SPACE  =  (  1 h?4,  (  <*30  ,101 ) 

// SYSLROD  00  DSN-ILGUSET(MAINI,UMT  =  SYSDA,01SP-(,PASS)« 

//  SPACE-ETRK, (10,10.?)  ) 

// 5 YSL  IN  DD  DSN=HK30)4.LKED1N4(LUT|,DISP=SHR 
II OJT  EXEC  PGM- •.LK£D.SYSlMCD,COND=  IB.LT  ,LKED) 

//FTOIFOOI  CO  UNIT  =  VI(J,DISP-(,PASS),SPAC£=ICYL,(),1)) 

//FT  0?F  001  DO  DSN-CUU, DISPMCJIO, DELETE  ) 

//FT03F001  00  DSN-LIF  IMP, DISPMULO,  DELI  TE  I 

//FT04F001  DD  OS  N-LOHPLS, 01  SP*(0LD, DELETE) 

//FTOSFDUl  DO  USN-U.F1FLD, D1SP-(0LD, DELETE  ) 

//FT06F001  DO  OSN *1 CPH 1  NT , D1 SP = ( M  U)  , PA S S) 

/  /  F  T  Of)  F  00 1  CD  USN-tCF  FIELD, 0ISP-(  (tv, DELETE) 

//FI09F001  DO  DSN-lLFPROO,Dl SP  = (OLD, CO) S) 

// FI10F901  UD  DSN-(.IRESP,01SP-(CL0»LD1S  ) 

// FT11F001  DD  DSN=CUA,DISP- (OLD.lDl  S) 

/If  1 1lf  001  CD  UNIT -V10,D)SP= I, PASS) , SPACE MCYL  ,( 1 ,1 )  ) 

//F  T  ?0F  00 1  CO  OS N - CLNL 1 N ,D ! SP - I GL D, 0 E L E TE  ) 

//  EXEC  PGM- IE  BGENER  »C  UN  0  =  E  VEN 

//SYS1N  00  DUMMY 
//SYSPRINT  DD  DUMMY 

//SYSJT1  DD  D5N-UPRINT  ,DJSP-(01D, DELETE  ) 

//SYSUT?  00  SYSOUT-tSYST, HOI 0-CHl CPRT 
//ODl  DO  DSN-llGOSET ,01 SP- (0L0,DELE  t  E ) 

//  EXEC  PGM-  TE  6CENE  R.CONO^E  VEN 

//SYSIN  00  DUMMY 
//SYSPRINT  DD  DUMMY 

//SYSJT1  00  DSN- CCPUNDAT A ,0)SP -(OLD, DELETE) 

//SYSUT?  00  SYSOUT -CPUN, HOLD -CHI PPL T 
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